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Abstract 

The average extensive air shower longitudinal development profile as a function of shower age 
in the energy range from 10^^ to lO^^eV is measured using data from the hybrid HiRes/MIA 
experiment. An angular bin signal based Cerenkov light component subtraction method and a 
shower maximum fitting method using a local parabolic function are used to correct and normalize 
the data. The Gaisser-Hillas and Greisen functions work equally well for describing the shower 
profile. The Gaussian function is a poor fit. A simple scale-free function is proposed and fits the 
data equally well. The best-fit parameters for the above functions are determined. 



I. INTRODUCTION 



The atmospheric nitrogen fluorescence hght technique plays an increasingly important 
role in extremely high energy cosmic ray observation. The success of the Fly's Eye and 
HiRes[0] experiments encourages the employment of this technique in future projects like 
the Auger experiment , the Telescope Array proposal^ and space based experiments like 
EUSO and OWL|p. In a recent report on the change of cosmic ray composition in 
the vicinity of 5x10^'' eV from HiRes/MIA joint dataQ, the authors have pointed out 
that the reconstruction of the shower longitudinal development depends upon the assumed 
functional form of the features of the shower development. How much do the characteristic 
parameters of the shower longitudinal development such as the shower maximum and its 
location, shower rise and shower decay constant depend on which function is employed? 
In the new technique of shower reconstruction developed for the HiRes monocular data 
0, the shower development function plays an even more fundamental role in the shower 
geometry reconstruction. In effect, the shower geometry is varied until the correct shower 
width predicted by the average shower function is obtained. The use of an appropriate 
shower development function is crucial in the determination of shower energy because this 
is based on the integral of this function. In practice, empirical functions based on data at 
lower energies or based on theoretical electromagnetic cascade calculation are used at the 
highest energies, e.g. above lO^^eV. However, none of these has been experimentally tested 
at these energies in the atmosphere. Values of parameters used in trial functions such as the 
A in the Gaisser-Hillas function and Lq (for a hadronic shower) in the Greisen function 
have never been measured in experiments at such a high energy. This is sufficient motivation 
to carefully test these well known functions and determine those parameters directly from 
experimental data. The HiRes/MIA joint experiment meets the necessary conditions for 
such a high resolution measurement of the shower longitudinal development. 

The HiRes prototype detector provides a degree by degree measurement of the intensity 
of the fluorescence light produced by the shower electrons along the shower axis. This 
light intensity, after corrections, is proportional to the shower development proflle. The 
angular bin signals, however, suffer fluctuations and are contaminated by the direct and 
scattered Cerenkov light component produced by the same electrons. The source of the 
fluctuation is sky noise and the statistical fluctuation due to flnite sampling of the light. A 



method to minimize the bin by bin fluctuation is to take averages of the angular bin signals 
over all the events. This will establish an almost perfect "average shower longitudinal 
development profile" which can be used to test the various trial functions. In order to 
do so, however, all showers must be "aligned" and "normalized" so that the showers can 
be compared with each other in spite of variations of position of shower maximum and 
shower energy. To carry out the alignment and normalization properly requires that the 
shower geometry must be well determined. In addition, the Cerenkov components of the 
signals must be subtracted. Fluorescence light is proportional to the shower size and the 
proportionality coefficient is geometry related. In the HiRes/MIA joint experiment, the 
shower geometry is well determined with the help of muon front timing information from 
MIA, the wide elevation angle range (up to 70°) of the HiRes prototype detector allows 
for a broad range of shower development to be detected, and the short shower-detector 
distance (about 3 km) minimizes uncertainties related to the atmospheric attenuation and 
scattering. We use the same data set selected for the study of the composition of primary 
cosmic rays||]. The criteria used to select this sample are well justified there based on a full 
detector simulation [p!0[] . 

In the following sections, we will describe the experiment and data set briefiy, explain 
how to minimize the contribution of Cerenkov light and how to subtract the Cerenkov light 
component in each degree bin, propose a way to normalize the individual showers, and 
present the average profile and compare it with well-known functions. 

II. HIRES/MIA JOINT EXPERIMENT AND DATA SET 

The HiRes prototype is situated at 112° W longitude and 40° N latitude and at a vertical 
atmospheric depth of 850 g/cw?. It is[|l| composed of 14 optical refiecting telescopes. They 
image the extensive air shower(EAS) as it progresses through the field of view of the tele- 
scope. Nitrogen fiuorescence light (300-400 nm) is emitted at an atmospheric depth X in 
proportion to the number of charged particles in the EAS at that depth. The electrons radi- 
ate Cerenkov light in the same band simultaneously. Because of the very foreward angular 
distribution of the Cerenkov light, only those events pointing toward the detector will give 
rise to signals dominated by Cerenkov light. However, the Cerenkov light can be scattered 
by molecular and aerosol particles into almost all pixels along the shower track no matter 



what the shower geometry is. Therefore, the signal recorded by the HiRes detector consists 
of the sum of fluorescence light, Cerenkov light and various kinds of noise. 

The MIA detector is situated about 3.3 km from HiRes and 150 m lower in elevation. 
It consists of 2500 m^ of active area distributed in 16 patches of 64 scintillation coun- 
ters and measures the EAS muon arrival times with a precision of 4 ns and records all hits 
occurring within 4 /is of the system trigger. Combining such well determined timing from 
MIA with the triggered HiRes pixel directions and HiRes timing information allows precise 
reconstruction of the shower geometry. This is essential in the determination of the atmo- 
spheric depth at which the shower passes the triggered HiRes pixels and the optical path 
length between the shower and the detector. 

HiRes/MIA coincident data were collected on clear moonless nights between Aug. 23, 
1993 and Aug. 24, 1996. The total coincident exposure time was 2878 hours corresponding 
to a duty cycle of 10.9%. 4034 coincident events were observed. For events passing a set 
of coincidence assurance cuts the shower trajectory, including arrival direction and core 
location for each event was obtained in an iterative procedure using the information from 



both HiRes and MIA ||T2[. The accuracy of the shower axis determination depends on the 
number of observed muons, the HiRes angular track length and the core distances from MIA 



and HiRes. 2491 events are reconstructed via this procedure. Monte Carlo(MC) studies fTO 



show that the median shower direction error is 0.85° with a median core location error of 
45 m. 

To insure data quality and maintain good resolution we require that for each event the 
depth of shower maximum Xm is visible within a minimum observed slant depth interval of 
250 g/cm^, that the track subtends at least 20° , that the accumulated gap between the fired 
pixels is less than 40% of the total gramage spanned by the shower, that the uncertainty 
in Xm, AXm, is less than 50 g/cm^, that the reduced fo^^ the profile fit not exceed 10 
and that the MIA to core distance Rp^iA ^^^^ than 2000 m. Additionally, we require 
the minimum viewing angle of the fired tube, 9^^ to be greater than 20° in order to delete 
events which possess large direct Cerenkov light components. These cuts, as summarized 
in TABLE 1, leave a sample of 488 events. The energy distribution peaks at 3x10^^ and 
the average Cerenkov light fraction is 15%. There is no event in this sample in which direct 
Cerenkov light contributes more than 1% of the total signal. 

We note that the shower development profile has been fit to the Gaisser-Hillas(G-H) 



function in order to carry out some of the cuts. However, the use of this function as well as 
the fit parameters derived from it is for the purpose of event selection only. Hereafter, only 
the raw data for these selected events is used in the analysis. 



III. CERENKOV LIGHT SUBTRACTION 

HiRes pixel signals are rearranged into a one dimensional series of one degree bins along 
the shower track by considering the light acceptance of the mirrors according to the geometric 
position of the shower. As mentioned above, these HiRes bin signals include both fluores- 
cence and Cerenkov light, although the total fraction of direct Cerenkov light is kept low by 
the imposition of the cuts described above. Bins which correspond to the deeper positions 
in the atmosphere tend to have more scattered Cerenkov light because the Cerenkov beam 
accumulates with increasing shower development and the molecular and aerosol scattering 
gets stronger as the shower penetrates deeper into the atmosphere. The raw signal versus 
depth, therefore, is not in itself a representation of the shower longitudinal development. 
We must subtract the scattered Cerenkov light from the observed signal. The resultant 
corrected signal composed only of the fluorescence light contribution is proportional to the 
shower size at the corresponding atmospheric depth except for the natural and electronic 
noise component. 

The Cerenkov light produced in an angular bin is calculated by multiplying the average 
Cerenkov light yield of a single electron with the shower size at the corresponding depth 
X. The average yield is taken over the energy of shower electrons above the Cerenkov 
threshold, Et, and over the emitting angle. The energy spectrum of shower electrons, the 
angular distribution of Cerenkov light and the pressure dependence of Ef and yield are 
carefully considered. Details can be found in the Ref. |]13[ and its references. 



The contribution of direct Cerenkov light is minimized by the data cuts. The contribution 
to the bin signal from scattered light from the accumulated Cerenkov beam is estimated 
as follows: we assume that the Cerenkov light beam in the angular bin is due to shower 
electrons in the previous adjacent angular bin. After taking into account the attenuation 
of the Cerenkov light beam traveling through the air between the bins, the contribution of 
scattered Cerenkov light from the beam into the detector is estimated based on the Rayleigh 
scattering theory and Mie scattering by aerosol. 



The observed signal in this angular bin can then be modified by subtracting the scattered 
Cerenkov light. The corrected signal can be converted into a shower size for the sake of 
computing the contribution to the Cerenkov light beam in the next adjacent angular bin. 
This recurrent procedure requires knowledge of the Cerenkov light beam component in the 
first angular bin in the field of the view of the detector. If the first signal bin were at the 
beginning of the shower development, the Cerenkov light component would be zero. The 
observed signal in this bin would be pure fluorescence which can be converted directly into 
a shower size. We assume that the flrst detected bin corresponds to an early enough stage 
in the shower development so that the zero-Cerenkov-light-component assumption is still 
approximately true. This subtraction scheme systematically underestimate the Cerenkov 
fraction in each subsequent bin because the first detected signal from the shower is not at 
the very beginning of the shower cascade. However, we will find that this systematic error 
is negligible and the proposed Cerenkov light subtraction is sufficient for our purposes. 



IV. THE NORMALIZATION OF AIR SHOWER DEVELOPMENT PROFILE 

As mentioned above, all the showers must be aligned in atmospheric depth and be nor- 
malized in shower size before one can take the average in bins over events. Most of our 
events follow a transition curve, namely the shower size increases rapidly as it develops in 
the atmosphere until reaches its maximum, Nm-, then it starts decreasing because shower 
electrons begin to lose more energy by ionization than by radiation of high energy gammas. 
The maximum size of a shower is approximately proportional to the energy of the primary 
particle. After subtracting the Cerenkov light components, the signals are converted into 
shower sizes by taking the geometry and atmosphere related light collecting efficiency and 
attenuation effect into account. The showers are normalized with respect to energy by their 
maximum size, i.e. N{X)/Nm, denoted as n{X). All the shower sizes are normalized to one 
as they reach their maxima. 

Air shower development fluctuates in atmospheric depth due to fluctuations in the 
hadronic multiparticle production (which depends on the nature of the primary nucleus) 
and electromagnetic processes. As a consequence, the position of the shower maximum, 
Xm, varies from event to event. One can express the shower longitudinal development as a 
function of the shower "age", defined as s = x+2X ' instead of depth, X. By using s, the 



shower longitudinal development is universally described as a rising phase from (initial 
position of the shower) to 1 (shower maximum) and the decay phase from 1 to 3 (infinite 
depth). The physical shower has an effective extent from to 2, however. 

After the Cerenkov light subtraction, the observed signal profiles of individual events 
span a range of atmospheric depth and show the shower maximum through their convex- 
upwards shape. The best way to find the position of shower maximum in the presence of 
statistical fluctuations is to carry out a local fit with a rather general function. We use a 
parabolic function to fit the data in order to determine both the location and the size of the 
shower maximum. The only constraint to this function is that it must be convex upwards. 
Here, "local fitting" means that the fitting procedure is only applied to those data points 
near the maximum to minimize the bias in the location of the maximum caused by using 
such a symmetric trial function. The method can be iterated to remove points far from the 
maximum. We find that the iteration affects the location of shower maximum by less than 
4% however and we do not use it in what follows. 

With shower maximum and size at maximum determined, the showers are normalized 
in size and aligned in depth. FIG 1 shows all 488 showers plotted with normalized shower 
size, n, versus shower age, s. In order to more easily see the change of the density of the 
dots in the scatter plot, we let each entry represents three measurements. 6069 entries in 
total out of 18230 bins signals are picked randomly and plotted. We note that there is a 
20% fluctuation in n near the shower maximum. The distribution of n at shower maximum 
is shown in FIG 2. The residuals in n, defined as the ratio between the deviation from the 
expectation value and the error associated with the individual bin signal, are well distributed 
as a Gaussian centered at with a width of 0.94. The main sources of this fluctuation are 
sky noise, statistical fluctuation, the error in shower geometry determination, the error in 
Cerenkov light subtraction and the error in determining shower maximum. The first three, 
and particularly the first one, contribute more than 80% to the total fluctuation in the 
signal. The error from Cerenkov subtraction contributes the least. 

The data is binned into age intervals, instead of the original angular bin. We take an 
average of the normalized shower size within each bin over all events to obtain an average 
shower development proflle. The curve is well determined between s=0.5 and 1.25. We 
demonstrate the effect of the Cerenkov light subtraction in FIG 3, by plotting the average 
shower transition curves before and after the subtraction in the same plot. 



V. TESTS OF TRIAL FUNCTIONS 



After calculating the average normalized shower longitudinal development profile, we fit 
the shape of the profile with several well known trial functions such as the Gaisser-HillasP|, 
Greisen0 and Gaussian forms. 

In order to describe the characteristics of the shower longitudinal development, i.e. the 
asymmetric rising-falling shape, Gaisser and Hillas introduced a function which depends on 
the essential characteristics an air shower: the initial point, Xq, the shower maximum, Nm, 
the shower maximum location, A^, and the shower decay length, A. Except for which is 
in number of shower electrons, all the parameters are in units of atmosphere depth, g/cm?. 
The function describes the change of the shower size with the atmospheric depth as 



N(X) = Njp^)~^e^, (1) 

where the atmospheric depth. A, is in g/cm?. 

Translating the depth A into age s and using the normalized shower size n = N/N^, 
eqn. (1) becomes 

where Tm = X^/A and Tq = Aq/A are the two remaining parameters. Parameter Tq is 
constrained to be less than „^'"'" Tm, where Smin is the lower limit of the data points (about 



rmn 



0.5 as shown in FIG. 3). The comparison of this function with the data is shown in FIG. 4. 
The two parameters are strongly correlated. 

The fit is poor at small and large s, particularly beyond s = 1.3. While a possible reason 
for the deviation of those 3~4 points may be poor statistics, a detailed systematic error 
analysis is described in the next section. 

The second trial function is the Greisen function which was suggested to describe 
the development of a pure electromagnetic shower. Depth is therefore expressed in the 
unit of radiation length, Lq = SQ.GQg/cm^. The function is defined by a single parameter, 
y = Xm/Lo, and has a form 

N{T) = ^L^e^(^-i'"^) , (3) 

where T = X/Lq refers to the atmospheric depth in radiation lengths in the air. Since the 
age s appears explicitly in the formula, it is straightforward to rewrite it as a function of 



either s or T. It is apparent after conversion to a function of s that the radiation length 
no longer appears explicitly in the formula. This feature of the Greisen function make it 
potentially useful for the description of hadronic shower with a single parameter y. However, 
the meaning of Lq need to be changed. In the case of electromagnetic cascade, Greisen shows 
that y = ln{Eo/Ec) fits a elongation law with a rate about 60g/cm^ where the Ec refers to 
the electron critical energy. We will see how the rate is changed to fit the case of hadronic 
shower below. The fit of the Greisen function to our measured average shower is shown in 
the dashed line in FIG. 4. 

In the Fly's Eye data analysis |14||, the poorer detector resolution dominated over the 



n{s) = expl-^ ( / '] ] }. (4) 



shape of the shower longitudinal development and a symmetric Gaussian in X could be 
used to represent the shower longitudinal development. Using ^ = cr/Xm as the single 
parameter, the Gaussian can be normalized by following the rule set in previous sections, 
i.e. 

The result is also plotted in FIG 4. It is a poor fit to the data with a = 26 per degree of 
freedom. Thanks to the significantly improved resolution, this experiment provides the first 
clear observation of an asymmetrical air shower development profile in X at 10^'' ~ lO^^eV. 

However, based on the result of this experiment, we observe that the average longitudinal 
development profile of air showers appears quite symmetrical as a function of shower age, 
s. A simple symmetrical Gaussian shaped function of age can be used to describe shower 
profiles. Since the function is centered at 1, it is governed by only one parameter, i.e. the 
width of the function, denoted as a. This fit is shown in FIG. 4 also. Hereafter, we call this 
function the "new Gaussian function" to distinguish it from the one that is symmetric in 
the depth variable. 

We plot the fit results in the form of deviations in FIG. 5. Quantitatively, the P^r 
degree of freedom for Gaisser-Hillas, Greisen and the "new Gaussian function" are 1.93, 
1.87 and 1.79 respectively. They all have essentially the same goodness of fit. More data 
beyond the region currently covered is required to further differentiate between them. From 
a practical standpoint however, the Gaisser-Hillas function contains too many parameters, 
especially Xq which is not directly measured in experiments and is consequently poorly 
defined. These parameters are strongly correlated. The Greisen function does not have this 



problem because it fits the data with a single parameter. Both functions need a scale to 
measure the shower development, however. The scale is given by A and Lq, respectively. 
The scale must be associated with the type of primary particle. This actually reduce the 
freedom of fitting in the reconstruction of real experimental events. One has to set the scale 
as a parameter in a complete fitting procedure. The "new function" proposed in this work 
is more practically useful in real shower reconstruction. As a function of depth, it has the 
form 



This function uses a single parameter, X^, to locate a shower in the atmosphere, and a 
single parameter, a, to indicate the shower width and it is scale-free. The last feature is 
very useful since the composition is a-priori unknown. 

VI. SYSTEMATIC ERROR ANALYSIS 

In order to investigate the reasons for deviation from the test functions at both early and 
late stages of shower development, we apply the above analysis to Monte Carlo simulated 
events. The generator used here is the full simulation code of the HiRes/MIA detector 
developed in the composition study It is driven by a shower generator based on the 
CORSIKA package. It produces the shower size using the Gaisser-Hillas function with full 
fluctuations in Xm, Nm, Xq and A. The fluctuation in muon density and muon arrival time 
on the ground is generated according to parameterizations of CORSIKA simulated events. 
All details of the detector are carefully considered including the sky noise, atmospheric 
attenuation, "ray-tracing", phototube response, electronics and triggering. Generated data 
are passed through the complete calibration and reconstruction procedure just like the real 
data. The details of the generator and the comparison with the data can be found in the Ref. 
||T0| . FIG. 6 shows that the overall quality of the fit is good. This means that the Cerenkov 
subtracting scheme, shower maximum finding and normalization procedure proposed in this 
paper work well. The shower longitudinal development function is reconstructed almost 
exactly except for the few points at the beginning and the end of the shower as is the case 
for the real data. The deviation turns out to be associated with the correction to tube signals 
which lie far from the shower-detector-plane. Details of this correction can be found in Ref. 




(5) 



. Briefly, a tube can be triggered even if it is not located exactly on the shower-detector 



plane since the shower has a finite lateral width, the mirror image on the focal plane has a 
finite spot size and the broadening effect due to the atmospheric scattering. The amount 
of light reaching the detector can be calculated by correcting the amount of measured light 
with a response function in the off-plane-angle. This function is determined by performing a 
careful "ray-tracing" procedure which folds all of the effects in. It turns out to be a Gaussian 
with a typical width of a few tenths of a degree depending on the shower distance and the 
shower age. On the other hand, we know that there is an uncertainty in the determination 
of the shower-detector plane direction of about 0.7 degrees pUj], This uncertainty causes a 
systematic overcorrection in the signals for off plane tubes. This becomes worse for a tube 
located far from the center of the phototube cluster. This effect can be verified by simply 
substituting the ideal MC input shower-detector plane for calculating the correction and 
comparing with the correction using the fitted shower-detector plane. This is shown in FIG. 
6. The effect on the beginning and end of the average profile is clear. 

One way to handle this systematic error would be to correct it based on this MC sim- 
ulation result. Another approach is to simply drop those affected points in the fitting 
procedures. We prefer the later for several reasons. First, it avoids introducing any model- 
ing dependence. Secondly, the points being dropped at age > 1.25, correspond to positions 
deeper than lASX^, or about 860g/cm^ in the atmosphere, where the uncertainty in atmo- 
spheric attenuation is greatest. As a part of the systematic uncertainty study, we varied the 
atmospheric parameters by one standard deviation from their most probable values in the 
shower reconstruction. We observed that there is no effect for most of the average profile 
and that this variation affected only the last few deepest points. This effect is smaller than 
that produced by the shower-detector plane error. Finally, in practice, dropping those points 
does not change our result because of the relatively large statistical errors associated with 
those points as can be seen in comparing the FIG. 4 with FIG. 5, where events associated 
with the last three points in FIG. 4 are dropped. 

We now discuss the systematic error caused by setting the Cerenkov light to be zero in 
the first angular bin in the field of view of the detector. Based on the previous results, 
we known that the Gaisser-Hillas function works well in describing the shower longitudinal 
development. We can use the Gaisser-Hillas function in a least x-square fitting procedure 
event by event. Since in this case we know the shape of the shower, the Cerenkov light 



component in all angular bins can be calculated. The fraction of Cerenkov light in each bin 
is compared to our original method in FIG. 7. It is clear that most of the Cerenkov light 
is subtracted with our bin signal based Cerenkov light subtracting scheme. Only 7.7% of 
the Cerenkov light is systematically underestimated. This is a minor effect in calculating 
the average transition curves as shown in the FIG. 3. The diamonds in the figure are 
obtained by carrying out the same averaging and normalization procedure as before but 
with the Cerenkov light fraction in each bin replaced by that determined using the Gaisser- 
Hillas function. The points are visibly lower which makes the profile appear slightly wider. 
However, the widening is not large enough to change our conclusion on the shape of the 
longitudinal development function. For instance, the values of fitted parameter Tq and 
are -2.57 and 8.6, respectively, which are consistent with the result listed below. 

The other systematic issue is the reliability of our maximum searching by using local 
fitting with a parabolic function. The uncertainty associated with this method is demon- 
strated by comparing the maxima with those obtained from the reconstruction procedure 
using the Gaisser-Hillas profile. FIG. 8 a) shows the comparison between the X^'s and FIG. 
8 b) shows the difference between the A^m's- The other way to illustrate the uncertainty is 
to compare the resolution functions in Xm and Nm from simulation for both methods. It 
turns out the resolution in for local parabolic fitting is 16% worse than the standard 
Gaisser-Hillas approach, and the resolution in is comparable in both cases. 



VII. RESULTS AND CONCLUSIONS 



In order to extract the best values for the parameters, T^, To, y and cr, in the Gaisser- 
Hillas, Greisen and the "new Gaussian function" , a well reconstructed data set is used. The 

parameters are determined to be 

T^ = 9.7±2.0, (6) 

To = -3.2 ±2.9, (7) 

y = 11.11 ±0.34, (8) 

a = 0.272 ±0.002. (9) 

A strong correlation between and Tq is found for the Gaisser-Hillas function. This 
correlation can be expressed as Tq ~ — O.IOT^ + 0.74T^ — 0.10 and is shown in FIG. 9. 



These results imply: 1) the strong correlation between parameters in the Gaisser-Hillas 
function means that the data can be fit with a simpler form. Alternatively one can fix one 
of those parameters and apply the correlation between them in the shower reconstruction; 
2) The parameter A in the Gaisser-Hillas function can be estimated to be 67 ± llg/cm^ 
by using the average Xmax = 651.9 ± 3.5g/cm^] 3) the parameter Tq, or Xq, is likely to 
be a negative number in any individual shower fit. It is thus difficult to interpret as the 
first interaction point of a shower. 4) the parameter y, defined as X^/Lq in the Greisen 
function, is close to T^, defined as Xm/ A. This implies that the best value of the parameter 
Lq, 58.7±2.0(7/cm^, is no longer the radiation length. This is to be expected since it is used 
to describe a hadronic shower; 5) a symmetrical function in age is a suitable representation 
of the shower longitudinal profile. This function is depth-scale invariant. 

In summary, the shape of the extensive air shower longitudinal development is investi- 
gated in the energy range from 10^^ to lO^^eV with the HiRes/MIA hybrid experiment. The 
profile is quite symmetrical as a function of the age of the shower. The "new Gaussian func- 
tion" @, the Greisen function (^) and the Gaisser-Hillas function describe the shower 
shape almost equally well, with the x^'s of the test 1.79, 1.87 and 1.93 respectively. This is 
the first direct measurement of the shower average longitudinal development profile at these 
energies and covering such a wide range of shower age. 
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Variables 


cuts 




Xi < Xjji < Xh 


Track Span 


Xh-Xi> 250 (g/cm^) 


irack Angular Length 


> 20° 


Total Gap in Gramage 


< QA{Xh - Xi) 




< 50 (g/cm^) 




<10/DOF 




> 20° 




<2 km 



TABLE L Event selecting criteria. Xh{Xi) refers to the depth corresponding to the highest (lowest) 
bin in the field of view of the detector. 
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FIG. 1: Aligned shower longitudinal development. Only one-third of the data is shown for clarity. 
Each point represents the fluorescence light signal in an one degree bin in an event. The events 
are normalized to 1 at individual shower maximum via a local parabolic fit. 
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FIG. 2: The fluctuation of normalized fluorescence light signals around shower maximum, 0.9 < 
s < 1.1. The dashed line represents a Gaussian fit, see the text for the details. 
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FIG. 3: The average shower transition curve. Sohd hne plus circles refers to the bin signal based 
recurrent Cerenkov-Hght-subtracting method, the dashed hne plus squares is the raw signal includ- 
ing the Cerenkov light. The dash-dotted line plus diamonds corresponds to a different Cerenkov 
light subtraction method (see the discussion on systematic error in Section V) . Lines are drown to 
guide the eye. 
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FIG. 4: Comparison between the data and test functions. Circles refer to data. The soHd hne 
refers to the G-H function, the dotted hne to the Greisen function, the dashed hne to the Gaussian 
function and the dash-dotted line to the newly proposed symmetrical Gaussian function of the 
shower age. 
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FIG. 5: Deviation of the data from the fits. The is per degree of freedom. 
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FIG. 6: The average longitudinal profile of MC generated showers. Open squares represent the 
results from the standard shower reconstruction; the open circles represent the same analysis but 
with the shower-detector plane direction used correcting the tube signal replaced with the ideal 
input parameters. 
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FIG. 7: Comparison of Cerenkov light fraction in each one degree bin. The x-axis is the recon- 
structed Cerenkov Ught fraction using G-H function and the y-axis is the fraction from the bin by 
bin recurrent Cerenkov-Hght-subtracting scheme. 
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FIG. 8: Comparison of shower maxima for reconstruction with the G-H function and with local 
parabolic fitting, a) for X^'s and b) for iVj^'s. 




FIG. 9: The correlation between the parameter and Tq in G-H function (|2|) 



